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We present a thermodynamic theory of the Kovacs effect based on the idea that the configurational 
degrees of freedom of a glass-forming material are driven out of equilibrium with the heat bath by 
irreversible thermal contraction and expansion. We assume that the slowly varying configurational 
subsystem, i.e. the part of the system that is described by inherent structures, is characterized by an 
effective temperature, and contains a volume-related internal variable. We examine mechanisms by 
which irreversible dynamics of the fast, kinetic- vibrational degrees of freedom can cause the entropy 
and the effective temperature of the configurational subsystem to increase during sufficiently rapid 
changes in the bath temperature. We then use this theory to interpret the numerical simulations 
by Mossa and Sciortino (MS), who observe the Kovacs effect in more detail than is feasible in 
laboratory experiments. Our analysis highlights two mechanisms for the equilibration of internal 
variables. In one of these, an internal variable first relaxes toward a state of quasi-equilibrium 
determined by the effective temperature, and then approaches true thermodynamic equilibrium as 
the effective temperature slowly relaxes toward the bath temperature. In the other mechanism, an 
internal variable directly equilibrates with the bath temperature on intermediate timescales, without 
equilibrating with the effective temperature at any stage. Both mechanisms appear to be essential 
for understanding the MS results. 
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I. INTRODUCTION 

The Kovacs effect reveals some of the most subtle and 
important nonequilibrium features of glassy dynamics. 
In particular, it provides detailed information about the 
ways in which glassy materials deform irreversibly and 
remember their histories of deformation Here, we 

develop a thermodynamic theory of the Kovacs effect, 
motivated in large part by the molecular-dynamics sim- 
ulations of Mossa and Sciortino (MS) [5]. 

In a Kovacs experiment, the volume of a glass- forming 
system is measured at fixed pressure as the temperature 
is varied. A sample is first quenched from a high tem- 
perature Th to a temperature low enough - i.e. near 
enough to the glass temperature Tg - that some inter- 
nal degrees of freedom fall out of equilibrium with the 
heat bath. The system is then aged at Ti, for times 
insufficient to reach thermal equilibrium, and finally is 
heated abruptly and held at a temperature Tf such that 

<Tf < Th- The crucial observation is that, in this 
last stage of the Kovacs protocol, the volume does not 
increase monotonically as a function of time, but goes 
through a maximum before decreasing slowly to its equi- 
librium value at Tf . The fact that the system can exist in 
two different states, on the upward and downward sides 
of the Kovacs volume peak, at the same temperature, 
pressure and volume, indicates that these are not states 
of thermal equilibrium. It is important to understand 
how to characterize them. 

The Kovacs effect originally was observed in polyvinyl 
acetate but since then it has been observed in many 
other glassy polymers, see for example @ for measure- 
ments in polystyrene. Qualitatively similar memory ef- 



fects have been observed in many other glassy systems 
such as colloidal glasses 0, ferroelectrics j^, fl0| - 
gelatin gels granular materials [l^l, superparamag- 
nets and superspin glasses [ll] . It also has been the sub- 
ject of various recent theoretical investigations [l^ - [20| . 

In this paper, we look at the Kovacs effect from the 
point of view of our recent attempts to develop a first- 
principles, statistical formulation of nonequilibrium ther- 
modynamics (2T| - [23j . Generally speaking, our goal is to 
reinterpret the analysis of Kovacs et al. [5] in terms of 
specific molecular processes. Our work differs from that 
of Nieuwenhuizen, Leuzzi, and coworkers [l^ [H, , for 
example, in that we start from a fundamental, statistical 
statement of the second law of thermodynamics and use 
it to derive equations of motion for relevant internal state 
variables as well as for an effective temperature. We dif- 
fer also from Berlin et al who have solved specific 
models and have shown how phenomena analogous to the 
Kovacs effect emerge in interesting ways. 

As in p2j . our starting point is the assumption that 
a glass-forming material consists of two weakly inter- 
acting subsystems. The configurational (C) subsystem 
is specified by the set of mechanically stable molecular 
positions, that is, the inherent structures [13, UHl- The 
kinetic-vibrational (K) subsytem is specified by all the 
other degrees of freedom - the kinetic energies, the dis- 
placements of the molecules from their stable positions, 
and, in the case of more complex molecules such as those 
relevant to the Kovacs effect, the internal degrees of free- 
dom of these molecules. The physical rationale for this 
separation is the distinction between the time scales for 
dynamic processes in the two subsystems. In situations 
where the system is driven by external forces, these two 



2 



subsystems can fall out of thermodynamic equilibrium 
with each other. The fast, K-subsystem remains in equi- 
librium with the heat bath at temperature T; the slower 
C-subsystem, at least transiently, has an effective tem- 
perature X (in energy units) that is different from ks T. 

We propose that any model of the Kovacs effect should 
include three essential ingredients. First, and most 
obviously, we need to specify the configurational (C- 
subsystem) degrees of freedom that fall out of thermal 
equilibrium as the whole system is quenched into the 
vicinity of its glass temperature. These degrees of free- 
dom must describe structural features that change via 
slow molecular rearrangements and thus do not keep 
up with more rapid variations of the bath temperature. 
Since we are interested in volume changes, the simplest 
choice of this internal variable is a population of vacancy- 
like defects; but other internal degrees of freedom that 
couple to the volume might serve our purposes equally 
well. 

Second, we need to ensure that the equations of motion 
for this out-of-equilibrium but statistically significant de- 
fect population are consistent with the laws of thermody- 
namics. We have argued in [22| that the natural way to 
do this is to use the effective temperature of the configu- 
rational subsystem, in direct analogy with Gibbsian sta- 
tistical mechanics, to determine the states of maximum 
probability through which the configurational degrees of 
freedom are moving. These configurational degrees of 
freedom have well defined energies Uc and entropies Sc] 
thus they have an effective temperature x = dUc/dSc, 
and their equations of motion must be based on their 
effective thermodynamics. 

Third, and least obviously, we need mechanisms by 
which variations of the ordinary temperature of the K- 
subsystem can produce changes in the effective temper- 
ature of the C-subsystem. This means that we must 
understand how, and under what circumstances, ordi- 
nary thermal expansion and contraction become irre- 
versible phenomena that can increase the entropy of 
the system as a whole. Our nonequilibrium thermo- 
dynamic formulation suggests that there are two dis- 
tinct thermo-viscoelastic mechanisms that can be rele- 
vant here. The first is a Kelvin- Voigt-type mechanism 
in which the K-subsystem exhibits a bulk viscosity aris- 
ing directly from fast molecular interactions. The second 
is a somewhat slower Maxwell-type mechanism involving 
volume-related internal degrees of freedom that equili- 
brate directly with the ordinary temperature T rather 
than the effective temperature x- (Our main reference 
for models of thermo-viscoelasticity is Maugin's book on 
The Thermomechanics of Nonlinear Irreversible Behav- 
iors, [261]). 

The work Mossa and Sciortino [5J offers a unique op- 
portunity to test the ideas described above. These au- 
thors performed molecular dynamics simulations of a Ko- 
vacs experiment using the Lewis and Wahnstrom model 
[27I I of ortho-terphenyl (OTP), in which the molecules are 
rigid isosceles triangles interacting via a Lennard- Jones 



potential. Their crucial result is that, when T is increased 
from to Tf, both the volume and the inherent struc- 
ture energy increase and go through maxima. This result 
tells us that the effective temperature and the vacancy 
population also increase and go through maxima, and 
that the thermal expansion driven by the change in the 
bath temperature is partially irreversible. More precisely, 
during thermal expansion or contraction at nonzero pres- 
sure, the system exchanges mechanical energy with its 
surroundings. Some of that energy is dissipated, produc- 
ing configurational entropy. 

One of the most remarkable results of Mossa and 
Sciortino [11 is shown in their Fig. 4. There they demon- 
strate that, after the volume and the inherent structure 
energy go through maxima, the system can be described 
by states of quasi-equilibrium fully characterized by the 
effective temperature. During the earlier stages, how- 
ever, this is clearly not the case. There, something else is 
happening that challenges our understanding of nonequi- 
librium thermodynamics. That "something else" is the 
central theme of the present investigation. 

The scheme of this paper is as follows. In Sec. HIl 
we describe our two-subsystem model and comment on 
its physical ingredients. The thermodynamic equations 
of motion for this model are derived in Sec. HID where 
we show how the two mechanisms of irreversible thermo- 
viscoelasticity emerge from our nonequilibrium statistical 
analysis. Section lTVl is devoted to identifying appropriate 
dimensionless variables and making first estimates of the 
parameters that appear in the scaled equations. 

In Sec. |Vl we compare the predictions of our theory 
with the numerical simulations of Mossa and Sciortino 
[^], and confirm that each of our three ingredients of 
a Kovacs model is, indeed, essential for understanding 
their data. During the reheating stage, as T rises quickly 
from Ti to T/, both the fast Kelvin- Voigt-type and the 
somewhat slower Maxwell-type mechanisms are needed 
as sources of configurational entropy. These sources in- 
crease the volume on short and intermediate timescales 
and drive an increase in the effective temperature x- The 
volume continues to rise as the vacancy population grows 
toward a quasi-equilibrium value determined by the in- 
creased X- After this quasi-equilibrium is established, at 
about the time that the volume reaches its Kovacs peak, 
the vacancies remain in equilibrium with x as the system 
slowly ages toward a final state with x — ksTf. In short, 
we recover the results of Mossa and Sciortino and fully 
agree with their interpretation of them. 



II. INGREDIENTS OF A TWO-SUBSYSTEM 
MODEL 

Denote the total, extensive, internal energy of the two- 
subsystem model by 

Utot = UciSc, Vel,N,) + Uk{Sk, Vel,Na). (2.1) 
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The configurational (C-subsystem) energy, Uc, is a func- 
tion of the configurational entropy Sc , an elastic volume 
Vei that is common to both subsystems, and an exten- 
sive number of vacancy-like defects Ny whose energies 
and excess volumes are e„ and Vy respectively. Similarly, 
the kinetic- vibrational (K-subsystem) energy, Uk, is a 
function of the kinetic-vibrational entropy Sk^ Vei, and 
the number Na of what we call "misalignment" defects 
with energies and excess volumes Va- 

Assume that a fixed number of molecules, say iVo, oc- 
cupies the elastic volume Vei , which does not include the 
excess volume of the defects. Then the total volume of 
the system is 



units) is 



Vel+NyVy+NaVa- 



(2.2) 



Our picture of the vacancy-like defects in the C- 
subsystem is a slight oversimplification but seems con- 
ceptually simple. In contrast, the misalignment defects 
require more discussion. The triangular geometry of an 
ortho-terphenyl molecule means that its volume and its 
energy depend on its orientation with respect to its neigh- 
bors. Thus, even in the absence of the vacancy-like de- 
fects that may characterize the slowly fluctuating con- 
figurational subsystem of OTP, there are local misalign- 
ment defects that couple to the volume and can partic- 
ipate in energetically irreversible processes. If the for- 
mation energies of these defects are not too large, and if 
the energy barriers that resist their transitions from one 
orientation to another are small enough, then these de- 
fects equilibrate quickly with the bath temperature and 
can legitimately be included in the kinetic-vibrational K- 
subsystem. For simplicity, we have assumed in Eqs. (|2.ip 
and (|2.2p that there is only one kind of misalignment de- 
fect. 

The assumption that the misalignment defects belong 
in the fast K-subsystem is not trivial. Our model is sim- 
ilar to one studied by Ilg and Barrat [ll] , who show that 
the equilibration rate for a similar class of dynamical 
inclusions in a driven glass former depends sensitively 
on the strengths of the thermal noise sources that acti- 
vate their transitions across internal barriers. We will 
show, however, that the high-temperature assumption 
works well for present purposes, and that these inter- 
nal, orientational degrees of freedom produce a model 
of bulk thermo-viscoelasticity that is consistent with the 
MS data. 

The temperature of the K-subsystem is 



Define the partial-pressure functions: 



and 



(dFc\ 

(dFK\ 



(2.4) 



(2.5) 



(2.6) 



where the free energies F are the Legendre transforms 
of the ?7's. Note, however, that we do not immediately 
identify pc+PK as the total applied pressure. We do, 
however, assume strictly linear elasticity by writing 



Fcix,Vei) Ac 



2Vg 



Vei-Vcix) +fc{x), (2.7) 



where Vo is a reference volume, Ac is a compression 
modulus, Vcix) is the relaxed C-subsystem volume, and 
fcix) is a volume- independent free-energy density. Sim- 
ilarly, 



Vo 



A 



K 



2Vi 



v,i-VK{e) +fKio). (2., 



The temperature dependent reference volumes Vc (x) and 
VK{d) are different from each other. The elastic energy 
of the C-subsystem has its minimum at a relatively small 
volume Vcix) J because that system is cohesive at zero 
pressure. In contrast, the energy of the K-subsystem 
decreases as the volume increases, because the kinetic 
energy is fixed and the vibrational modes become softer 
as the spacing between the molecules increases. Thus, 
strictly speaking, Vk{0) must actually be defined at a 
positive reference pressure; but, since we are considering 
only linear elasticity, there is no need to be specific about 
this definition. In any case, the partial pressures defined 
in Eqs. (|2.5I) and (|2.6p are different from each other and 
most likely have opposite signs, i.e. pc <0<Pk- This is 
indeed the case in the simulations of MS [al . 



III. THERMODYNAMIC EQUATIONS OF 
MOTION 



First and Second Laws 



K 



\dSi. 



(2.3) 



We assume that the K-subsystem, in addition to having 
its own internal dynamics, plays the role of a thermal 
reservoir, so that 9 is the temperature that is being con- 
trolled as a function of time during a Kovacs experiment. 
The effective temperature of the C-subsystem (in energy 



The first law of thermodynamics is 
- pVtot = Utot, 



(3.1) 



where p is the applied pressure. With Eqs. (|2.2p . (|2.3p . 
and (EH), Eq.dSU) becomes 



xSc 



pvy 



dUc 



Ny+ 9Sk 



4 



pVa 



[ OUk 



Na 



+ [p-Pc{x,Vei)~PK{d,V,i)\V,i^O. (3.2) 
The second law is 

Stot — Sc + Sk> 0. (3.3) 



Using Eq. (|3.2l) to eliminate Sc, we write Ea. (l3.3p in the 
form: 



PVa 



pvo + 



fdU^ 



K 



\dNa 
(dUc\ 



Sk,V^i 



Na 



(3.4) 



P~Pc{X,Vel) -PK{0,Vel) 



Vel 



{0 ~x)Sk> 0. 



Following the procedure described in [2l|, [24I, we rec- 
ognize that Eq. (|3.4p consists of four separate inequalities 
associated with the independently variable quantities iVa, 
Ny, Vel, and Sk, and therefore we must satisfy four sep- 
arate inequalities. In the next paragraphs, we look at 
these in reverse order of their appearance here. 



B. Aging Rate 



The inequality proportional to Sk is satisfied by writ- 



ing 



es 



K 



-Aix,0) (l-f) 



(3.5) 



where A{x, 9) is a non-negative thermal transport coeffi- 
cient. 9 Sk is the rate at which heat is flowing from the 
C-subsystem to the K-subsystem. Since we assume that 
the coupling between the C and K-subsystems is weak, 
we expect A{x, 9) to be small. This is the term that con- 
trols the rate at which the system ages in the absence of 
external driving. 



C. Kelvin- Voigt- Type Thermo- Viscoelasticity 



The inequality in Eq. (j3.4p is satisfied by writing an equa- 
tion of motion for V^i '■ 



To Vel 



'7 



Vel 



Ve7{x,9,p) 



(3.8) 



where tq is a molecular time scale and tq /j is a bulk vis- 
cosity. According to Maugin [2^, this is a Kelvin- Voigt- 
type thermo-viscoelasticity. When rewritten in terms of 
the pressures, Eq. p.Sp says that the driving force p is 
equal to an elastic term, pc+Pk, plus a viscous force 
proportional to Vei ■ 

To see in more detail what is happening here, interpret 
the first-law in Eq. p.2p as an equation for x Sc, and note 
that the contribution to the configurational heating rate 
from the term proportional to Vei is 



P-PC{X, Vel)-PK{9, Vel) 



Vel 



Vo 



[Vel 



-v:i\x,^ 



\p) 



(3.9) 

Ordinarily, this term is negligible. In the absence of a 
slow, internal, dissipative mechanism, the viscosity ro/7 
is microscopically small. Suppose that some quantity on 
the right-hand side of Eq. (|3.8p . perhaps p or 9, is var- 
ied at an experimentally feasible rate, say v/tq <^^/tq. 
By dimensional analysis of Ea. (l3.8p . we find that v/j^ 
{Vei — V^^)/Vo. If we then integrate the right-hand side 
of Eq. (|3.9p over a time of the order of , we find that 
the change in the total heat energy is of the order of 
Wqv/^, which vanishes when 0. In this limit, 

the system becomes thermodynamically reversible. The 
solution of Ea. p.8p is accurately p — Pk+Pc, which is 
the usual thermodynamic identity. In other words, we 
have recovered a special example of the general rule that 
equilibrium thermodynamics is valid when systems are 
driven quasistatically. 

But the Kovacs effect is an exception to this rule. The 
original Kovacs observations were made with a glassy 
polymer, where the internal timescales ro/7 are long, so 
that it is possible to change temperatures and pressures 
relatively rapidly. It is also quite easy to do this in molec- 
ular dynamics simulations, which is what happens in the 
MS computations. We will see in Sec. |V]that the Kelvin- 
Voigt-type dissipation is an important driving force for 
the Kovacs effect. 



Next, consider the part of the inequality proportional 
to Vel- Use the definitions of the partial pressures in Eqs. 
p.Sp and (12. 6p . plus the elastic free energies in Eqs. (|2.7p 
and (|2.8p . to write 



p-pc{X,Vel)-pK{9,Vel) = ^ [Vel-V:i\x,9,p) 



where A = + Xk, and 



(3.6) 



Verix,9,p) 



XKVK{9) + XcVcix)-pVo . (3.7) 



D. Maxwell- Type Thermo- Viscoelasticity 

We turn finally to the terms proportional to Na and 
Ny in Eq. (|3.4p . These terms produce a Maxwell- type 
thermo-viscoelasticity, according Maugin [2^. To see 
this, look at the time derivative of the expression for 
the total volume in Ea. (|2.2p . If there is no Kelvin- Voigt- 
type viscous pressure proportional to Vei , then the total 
deformation rate Vtot is the sum of an elastic term Vi — 
— Vbp/A, and a viscoelastic term Va Na + v^N^. Our 
problem reduces to finding an expression for Va Na+Vy Ny 
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that will serve as a constitutive relation for the viscoelas- 
tic (viscoplastic) part of the deformation rate. To solve 
this problem, we follow steps outlined in [2ll |. 

Start with the K-subsystem. Assume that the entropy 
and energy of this system consist of separate, additive 
contributions, first from the defects and, second, from all 
the other degrees of freedom: 

SK{UK,V,uNa) = S^iNa) + SiiUuVei); (3.10) 



UK{SK,VeuNa) ^ Naea + Ui{SuVel)- 



(3.11) 



where, for simple defects without internal structure of 
their own, 

Sq{N) = TVo IniVo - N\nN -{No-N) \tl{Nq ~ N), 

(3.12) 

and A^o, as defined earlier, is the total number of 
molecules. Then, 



UKiSK,Vel,N,)=Naea + Ui 

SO that 



'if 



K 



dNa 



Sk,V^i 



dN„, 



(3.13) 



(3.14) 



The inequality associated with the Na term in Eq. (13. 4p 
has the form of a Clausius-Duhcm relation ^] , enforcing 
non-negative entropy production: 



pva + 



fdU, 



K 



\dNa 



dG 



K 



(3.15) 



where 

GK{P.V.Na)^haNa~BS^{Na)\ (3.16) 

and ha = Sa +pva- We satisfy Eq. p.lSp by writing 

To Na - Tk [K'iO,?) - Na] , (3.17) 

where Tx is a dimensionless rate factor and tq is the 
same molecular time scale that we introduced in Eg. p.Sp . 
Na'^{0,p) is the equilibrium value of Na determined by 



fdG 



K 



V dNa / 
which means that 



= at Na^K'io,p); 



No 



(3.18) 



(3.19) 



The same analysis pertains to the vacancy-like defects 
in the C-subsystem. The equation of motion for Ny is 
the same as Eg. (13. 171) . but with v-u and replacing Va 
and Ba, with x instead of 6, and with a new rate factor 



ToNy=Tc [N^^ix,p)- Ny] 



(3.20) 



We assume that the entropy associated with Ny is the 
same function So{Ny) that we introduced in Ea. p.l2p : 
therefore 



N^'ix,P) 



No 



ehjx + 1' 



(3.21) 



where /i„ = Cy + pVy. We note that the idea that 
an internal variable can be transiently driven out of 
quasi-equilibrium with the effective temperature, as in 
Eg. p. 201) . was introduced earlier in [2^, where it was 
used to describe the internal dynamics of deforming, sim- 
ulated, amorphous silicon. 



E. Equation of Motion for the Effective 
Temperature x 



Putting these pieces together, we rewrite Eg. p. 21) as 
an expression for the heat flow into the C-subsystcm: 



X Sc- 
7A 



-X- 



dSoiNy) 



dNy 



Ny^ 



ha - 



dSoiNa) 



dNa 



Na 



+ ^—[V,i-V:,%x.0.p)\ +A{x,6) (l-|) 



(3.22) 



The first three of these terms are non-negative rates of 
configurational heat production; the last term is the (or- 
dinarily negative) rate at which heat flows from the K- 
subsytem into the C-subsystem. 

To convert this result into an equation of motion for 
X, write 



xSc 



dSc . _^ dSc 

X -TT- X + X TTT^ N, 



dx 



dNy 



dSc 

X -^77- Vel 

dVel 



= + ,3.3, 

The second term exactly cancels the term proportional 
to d So /dNy on the right-hand side of Eg. p. 221) . which 
becomes 



C'ff X=-hyNy 



dpc 

dx 



ha 



dSoiNa) 



dNa 



Na 



+ ^ [V,i^V:,''ix,0,p)]\Aix,9) (l-|) . (3.24) 



IV. SCALING AND APPROXIMATIONS 

We have arrived at a complex set of equations with 
many variables and parameters. Before using these equa- 
tions for data analysis, we rewrite them in terms of di- 
mensionless quantities and, where possible, identify phys- 
ically motivated estimates for some of the parameters. 

To start, rescale the time so that tq = 1. Then define 
the defect densities: 



Na, 

No' 



(4.1) 
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and the volume fractions: 
, Vtot , 

9tot = -Jl—i Vel = 
1^0 



(4.2) 



We also need to express the defect volumes in units of 
the volume per molecule: 



No 



Vy-, (t>a 



No 

Vo 



Therefore 



and 



</>tot = (l)el + (l)ana+ 4>v riy] 



(pel = -1 [(Pel 

Measure x in units of the vacancy enthalpy: 

~ = 2L- 

hy 



(4.3) 

(4.4) 
(4.5) 
(4.6) 



thus, 



1 



;1/X + 1 



(e-'/*-n,), (4.7) 



where the last approximation is valid in the low-density 
limit, 1- For comparison with experimental data, it is 
convenient to express 9 in units of absolute temperature 
T, so that 



K 



1 



(4.8) 



where ks Ta = ha. 

For simplicity, assume that the relaxed volume of the 
C-subsystem, Vcix)^ is independent of x, and write 



0:«(T,p)-0o+0i(T), 



(4.9) 



where (j)o — 1— p/A, and (pi{T) describes the ordinary 
thermal expansion and contraction that drive the Kovacs 
experiment. In general, (j)i{T) is a nonlinear function 
over the range of temperature jumps used by MS, and 
we will need to use their data to evaluate it. 

The equation of motion for x, i-e. Ea. p.24p . now reads 



T , 
IH In , 



+ 76 (pel- (pJiT,p) 



rA(x,r) ( _ _ ) . (4.10) 



Lennard- Jones potential, we estimate the molecular vi- 
bration period to be about 1.5 picoseconds. Therefore, it 
is convenient to choose our unit of time to be tq = 1 
ps. Similarly, from the characteristic energy scale of 
this potential, we estimate that hy ~ 0.1 eV, so that 
Ty ^ 1300 K. The misalignment defects must have sub- 
stantially smaller formation energies. If we guess that 
the difference is roughly a factor of ten, then 130 K; 
so that these defects are far from being frozen out at the 
lowest temperatures used by MS, i.e. at = 150 K. This 
same estimate of hy, combined with Vq/Nq ^ 0.4 nm'^, 
and an estimate for the bulk modulus A ^ 4 GPa, tells 
us that b ~ 50, which means that the heating rate asso- 
ciated with the Kelvin- Voigt- type term in Eq. (|4.10|) may 
be substantial. 

These estimates have interesting implications for our 
choices of the rate factors. Clearly, with Ta~130if, the 
transition rate for the misalignment defects is not appre- 
ciably limited by an activation barrier; so Tk must be 
only moderately slower than the molecular rate 7, which 
by definition cannot be significantly different from unity. 
Our first guess is that is in the range 10~^ to 10~^. 
On the other hand, Tc should contain an effective ther- 
mal activation factor of the form exp(— Ac/x), where 
Ac is the excess barrier, in units of hy, that the system 
must surmount in either creating or annihilating a va- 
cancy. If we assume that the system is fully equilibrated 
at T = Th = 400 K, then the initial value of x is equal 
to ksTfi/hy '^0.53. Assuming that Ac is not too much 
smaller than unity, we conclude that Tc may be smaller 
than Tk by a factor of ten or so. 

The more interesting rate factor is Ta, which, un- 
like Tc in Ea. (|4.7[) . is not the prefactor in a creation- 
rate formula that already contains an activation factor 
exp(— 1/x). In other circumstances, such as a calcula- 
tion of the a relaxation rate or the shear viscosity in 
a glass forming material that is moving so slowly that 
X~kBT, the analog of Ta would be a super-Arrhenius 
function of T. Here, although we are talking about the 
slow aging part of a Kovacs experiment, we are looking at 
the early transient stage where x is still somewhat bigger 
than T/Ty in Ea. (|4.10p . Accordingly, we assume that this 
rate at which the two weakly coupled subsystems equili- 
brate with each other is limited by a substantial energy 
barrier, say A a > 1; and we write 



r^(x,T)-r^°)(T)e-^-/*. 



(4.11) 



Here, d'^f = C'^^ /No, ke Ty ^ hy, p ^ ha/hy = Ta/Ty, 
6= (A Vo/hy No) and r^(x, T)^A{x, 6) /{No kg T). 

We can make rough estimates for some of these pa- 
rameters using known properties of the OTP model sim- 
ulated by MS. For example, taking parameters from their 



If the super-Arrhenius analogy is valid, F^^ (T) will be a 
rapidly varying function of T that becomes vanishingly 
small below Tg. At constant T, while x>T /Ty, Eq. fiTT)) 
implies that the aging rate slows exponentially as x de- 
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V. COMPARISONS WITH THE DATA OF 
MOSSA AND SCIORTINO 

In their molecular dynamics experiments, MS are able 
to resolve dynamic variations in their model of ortho- 
terphenyl on time scales as small as tens of picoseconds. 
In addition to observing volume changes, they can ob- 
serve the energy, the pressure, and the shape factor (a 
measure of the width of the energy basins) of the inher- 
ent structures throughout the Kovacs protocol. Thus, 
they probe the Kovacs phenomena to a depth that seems 
impossible for laboratory experiments. 

There are, however, compromises that must be made in 
such a procedure. MS simulate a system of only 343 OTP 
molecules. Although they average their results over hun- 
dreds of initial configurations, it is hard to rule out effects 
of numerical noise, especially in the low-temperature ag- 
ing calculations that must be dominated by very rare 
events. 

Moreover, to control temperature and pressure, MS 
use a thermostat and a barostat with a time constant of 
20 ps; and they state that their systems are too far out 
of equilibrium for the data to be meaningful on shorter 
time scales following the initial quench or the reheating 
step. In spite of these uncertainties, we decided to try 
to model the complete Kovacs data reported by MS. We 
computed the values of the volume, the effective temper- 
ature, and the defect densities at the end of the aging 
stage, and used these values as the initial conditions for 
the reheating stage. Note that the interesting features of 
the Kovacs effect are very small; i.e. the fractional vol- 
ume change associated with changes in x near the Ko- 
vacs peak is only of the order of 10"'^. Therefore, in some 
places, we have adjusted the values of our parameters to 
three or more significant figures in order to make quan- 
titative comparisons with the MS data. 

Our data fitting procedure for the results presented 
here started by choosing 

(jQ{T = 400 if ) = 00 = 1 - p/A = 0.9947, (5.1) 

which was based on the estimate p = 16 MPa and A — 
3 GPa. We then estimated r„ = 1300 K, Ta = 130 K, u„ = 
0.07 nm^, and = 0.007 nm^. Assuming full equilibrium 
at Th — 400 if, X — Tfi/Ty, we computed Uy and Ua at 
that temperature. MS report that their total volume per 
molecule at 400 if is 0.378 nm'^. These numbers uniquely 
determine Vq /Nq = 0.374 nm^ . 

From here on, we chose parameters in accord with our 
rough estimates at the end of Sec. IIVI and refined these 
estimates to improve the agreement with the data. In 
addition to those cited in the last paragraph, the fol- 
lowing numbers were used throughout the calculations: 
b = 30, Aa = 3.5, Tk = 10-l■^^ and Pc = 10" ^-^s^ 
Our best-fit values for the elastic volume fractions at 
T=150if and 280 if were (T = 150 if ) = 0.9037 and 
(T = 280 if) = 0.9323. 

In Figs. [Hand [21 we show theory and data for the MS 
instantaneous quench from Th = 400 if directly to Tf = 
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|°9io(y 

FIG. 1: Aging at Ti — 280 if after an instantaneous quench 
from Th = 400 if. Upper panel: time evolution of the theoret- 
ical volume per molecule Vo/No (solid line) compared to the 
simulation data (open circles) extracted from MS Fig. 1(c) 
3. The parameters used for the theoretical curve can be 
found in the text. Lower panel: the corresponding reduced 
effective temperature x^x/^v- 
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FIG. 2: Ua, riv and (j>ei corresponding to Fig. [1] 



280 if and subsequent aging as functions of the time after 
quench te- (AH times are stated in picoseconds.) The 

theoretical parameters were 7 = 10^ ^'^^ and P^"* = lO"'^'*. 
We need a Kelvin- Voigt-type viscosity with a small 7 
because the elastic part of the volume initially relaxes on 
a time scale of the order of 7" ^ ^ 15 ps. The resulting 
dissipation drives a rapid increase in and then both 
Ua and Uy participate in the change of the total volume 
on time scales determined, respectively, by Tk and Pp. 

Next consider the quench from = 400 if to = 
150 if and subsequent aging at the latter temperature. 
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FIG. 3: The same as Fig. [U but for Te = 150 K, after a 
smooth quench from Th = 400 if according to EQ. (|5.2p . See 
text for the parameters used. The simulation data extracted 
from MS Fig. 1(c) H. 



In this case, MS have told us that they used a smooth, 
thermostaticahy controlled decrease in the temperature, 
and started to measure the volume at about 30 ps after 
the quench was started [sO] . Accordingly, we have shifted 
our time scale by 30 ps; and we have modeled the initial 
temperature dependence by writing 



T{te) = T, + in - Tt) exp ( 

Tth 



(5.2) 



with Tj;j =4ps. To use this equation for values of T be- 
tween Th and T^, we have made a linear interpolation of 
(1)11 (T) between the values given above for 4f^{Th) and 
(jffiTi.). We used 7 = 1 and = IQ-o ^. The results 
are shown in Figs. [3] and 21 along with the MS data 
for the volume. With the more gradual quench and the 
larger value of 7, the Kelvin- Voigt-type effect is less pro- 
nounced but still present. The important feature here 
is the very slow aging at long times, associated with 
the smaller value of F^ and the controlling effect of the 
effective-temperature activation barrier A^i. 

The principal Kovacs effect occurs during and after 
reheating from T^^ 150 if to T/^ 280 if, starting at the 
end of the aging period t = te- Because MS increased T 
smoothly during reheating, we have used 



Tit)=Tf-iTf-T,)exp l- 



t - 

Tth 



(5.3) 



again with Tth = 4 ps, and with a linear interpolation 
between the lower and upper values of (j/J{T). For this 

stage, we have used the same values of 7 and fJ^'' that we 
used for aging at 280 K. The results for a waiting time of 
te — 25 ns (equivalent to the nominal MS value of 25 ns on 
our shifted time scale) are shown in Figs. [5]and[6l along 
with the MS data for the volume. Again, we find that we 
need all three irreversible mechanisms to understand the 
observed behavior. The effective temperature increases 
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FIG. 4: 



and (j>ei corresponding to Fig. |31 



quickly, due to the combination of both Kelvin- Voigt- 
type and Maxwell-type mechanisms. The density of K- 
subsystem defects, n^, rises toward equilibrium with Tf 
at a rate Tk', then the vacancy density Uy rises toward 
equilibrium with % at a rate Fc; and finally Uy, now at 
its quasi-equilibrium value as a function of %, decreases 
as X decreases slowly toward ksTf. 

Figure [7] shows the Kovacs peak from Fig. [S] and the 
comparable peak for a shorter waiting time, = 1 ns, 
along with the MS data for both cases. The agreement 
seems excellent in view of the fact that we computed 
the second curve only after having determined all of the 
parameters from the preceding calculations. 

Finally, in Fig. |51 we show the inherent-structure en- 
ergy ejs for the reheating stage shown in Fig. O To fit 
the MS data, we have used 



e/s 



(pel ~ (PlU^i) 



naks Ta + iiv ks Ty 

(5.4) 



with e/s = -84.32 kJ/mol and e'jg = 0.08. Note 
that our inherent-structure energy e/5 contains not only 
Uc{Sc,Vei, Ny), but also a contribution from the mis- 
alignment defects Na that we earlier argued belong to 
Uk for thermodynamic reasons. The conventional, static 
definition of the inherent-structure energy requires this 
interpretation of e/s- We use the canonical variables T 
and X in Ea. (|5.4p instead of the entropies that we used 
in the micro-canonical formulation in Eq. (j2.ip . Inter- 
estingly, the position of the peak in ejs{t) seems to be 
determined most strongly by the time dependence of the 
configurational vacancy term, i.e. HyksTy, as opposed 
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FIG. 5: The same as Fig. [U but for reheating from Ti = 150 K 
to Tf — 280 K according to Ea. (|5.3p . and after an aging time 
logio(ie) = 4.4 (equivalent to the nominal MS value of 25 ns 
on our shifted time scale). See text for the parameters used. 
The data in the upper panel are extracted from MS Fig. 2 
. The Kovacs peak appears to be small because the vertical 
scale includes the initial thermal expansion. See Fig. [7] for a 
closer look at the peak. 
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FIG. 6: and 4>ei corresponding to Fig. O 



to being determined predominantly by other configura- 
tional degrees of freedom and thus more directly related 
to the time dependence of x- 



VI. CONCLUDING REMARKS 

In their own concluding remarks, Mossa and Sciortino 
summarize their results by saying that, instead of 
moving a system along a sequence of quasi-equilibrium 



FIG. 7: The same as the upper panel of Fig. O but after an 
aging time logj^o(fe) = 3 (equivalent to the nominal MS value 
of 1 ns on our shifted time scale). The upper panel of Fig. [5] 
is copied here for comparison. See legend for more details. 



configurations, their "aging dynamics propagates the sys- 
tem through a sequence of configurations never explored 
in equilibrium, and it becomes impossible to associate the 
aging system to a corresponding liquid configuration." 
They go on to ask whether "a thermodynamic descrip- 
tion can be recovered [by] decomposing the aging system 
in a collection of substates, each of them associated with 
a different fictive T ... or if the glass ... is trapped in 
some highly stressed configuration which can never be 
associated with a liquid state." 

We seem to be arriving at a related but different inter- 
pretation. By focusing on internal state variables - in this 
case, the density of different kinds of defects in both the 
configurational and kinetic-vibrational subsystems - in 
addition to the effective temperature, we naturally gen- 
erate states in which the system as a whole departs from 
both ordinary thermal equilibrium with the bath tem- 
perature and from quasi-equilibrium with the effective 
temperature. Our technique for making this calculation 
is the one we described in [2l|, . We think that this 
technique goes at least part of the way toward answer- 
ing the questions posed by Mossa and Sciortino, but we 
recognize that it encounters conceptual problems that 
eventually must be addressed. 

The most obvious such problem, in our opinion, is the 
one that we found when deciding to include the misalign- 
ment defects among the fast degrees of freedom in the 
kinetic-vibrational subsystem. The results of Ilg and Bar- 
rat f2§\ imply that some such internal variables may be 
neither completely fast nor completely slow but, rather, 
their dynamics might be activated by an intermediate 
temperature or noise strength. This possibility might 
be loosely related to the MS conjecture about "different 
fictive [temperatures];" but the conjectures are intrinsi- 
cally different from each other. Neither we nor Ilg and 
Barrat are contemplating more than one effective tem- 
perature. So far as we can tell, no complication of either 
kind is needed for understanding the Kovacs effect as ob- 
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FIG. 8: The theoretical inherent structure energy e/s of 
Eq. (|5.4|) (solid line) compared to the simulation data (open 
circles) extracted from MS Fig. 3(b). 



served by MS, nor do we seem to need it for shear flow in 
amorphous systems [1^ or in polycrystals [sij . Neverthe- 
less, we appear to be encountering some of the deepest 
and most important open questions in nonequilibrium 
physics. 
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